Quantization of the chiral soliton in medium 
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Chiral solitons coupled with quarks in medium are studied based on the Wigner-Seitz approxima- 
tion. The chiral quark soliton model is used to obtain the classical soliton solutions. To investigate 
nucleon and A in matter, the semi-classical quantization is performed by the cranking method. The 
saturation for nucleon matter and A matter are observed. 
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I. INTRODUCTION 

The study of dense nuclear matter with the internal 
nucleon structure is old but still a challenging subject. 
Especially, the approach of the topological soliton model 
seems interesting, because it is believed as a low energy 
effective model in the large iVc-liniit of QCD. It was first 
applied for nuclear matter system in 80 's by using the 
skyrmion centered cubic (CC) crystal by Klebanov 0. 
This configuration was studied further by Wiist, Brown 
and Jackson to estimate the baryon density and discuss 
the phase transition between nuclear matter and quark 
matter Goldhabor and Manton found a new con- 
figuration, body centered cubic (BCC) of half-skyrmions 
in a higher density regime Q- The face centered cubic 
(FCC) and BCC lattice were also studied by Castillejo et 
al. 4] and the phase transitions between those configura- 
tions were investigated by Kugier and Shtrikman Q . Re- 
cently, the idea of using crystallized skyrmions to study 
nuclear matter was revived by Park, Min, Rho and Vento 
with the introduction of the Atiyah-Manton multi-soliton 
ansatz in a unit cell Q. 

Incorporating quark degrees of freedom into each soli- 
ton makes the prediction more realistic. Achtzehnter, 
Scheid and Wilets investigated the Friedberg-Lee soli- 
ton bag model with a simple cubic lattice Due to 
the periodicity of the background potential, the solution 
of the Dirac equation has the form of the Bloch waves, 
fpkir) — e^''"^ 4>k{r) where (j)k satisfies the same periodic 
boundary condition as the background potential. They 
performed the calculation for only one direction of the 
crystal momenta k = kez and assumed the spherically 
symmetric energy surface. The Bloch condition is, how- 
ever, anisotropic for the nonzero k and the results should 
be highly dependent on the approximation. The analy- 
sis of the crystal soliton model with quarks based on the 
Wigner-Seitz approximation has been already done. In 
this ansatz, a single soliton is placed on the center of a 
spherical unit cell. Then the lowest energy level ("bot- 
tom" of the band) for the valence quarks becomes s-state. 
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The appropriate boundary conditions at the cell bound- 
ary should be imposed on the quark wave functions as 
well as the chiral fields. This simple treatment sheds light 
on the nucleon structure in nuclear medium. Soliton mat- 
ter within this approximation have been extensively stud- 
ied by using various nucleon models such as the the chi- 
ral quark-meson t ype model lEJM ITol ITlj , Friedberg-Lee 
soliton bag model [llllllIlQrthe Skyrme model [H. 
The non-zero dispersion of the lowest band and the 
quark-meson coupling [l^ were also examined within this 
approximation. 

The chiral quark soliton model (CQSM) can be inter- 
preted as the soliton bag model including not only va- 
lence quarks but also the vacuum sea quark polarization 
effects explicitly ^16,, 17, 18]. The model provides cor- 
rect observables of a nucleon such as mass, electromag- 
netic value, spin carried by quarks, parton distributions 
and octet, decuplet SU{3) baryon spectra 19]. Remark- 
ably this model predicted the exotic quark bound state, 
pentaquark j20| which was successfully observed in 
experiments [2lj] . 

Amore and De Pace studied soliton matter in the 
CQSM using the Wigner-Seit z ap proximation and ob- 
served the nuclear saturation [23]. They examined the 
soliton solutions with three different boundary condi- 
tions imposed on the quark wave function. However the 
obtained saturation density was higher than the experi- 
mental value and they concluded that such discrepancy is 
originated in the approximate treatment of the sea quark 
contribution . Thus we treat the vacuum polarization 
exactly in the manner originally proposed by Kahana and 
Ripka .24] and semi-classically quantize the chiral soliton 
by the cranking method to see those effects on the mat- 
ter solution. At present, soliton matter has been studied 
only at the classical energy level. In order to study the 
property of nucleon or A in medium, the spin and isospin 
of each of the soliton must be quantized. We hence per- 
form the rotational collective quantization by the crank- 
ing formula and observe the saturation of nuclear and A 
matter. As shown in Sec. V, we obtained very shallow 
saturation. 

Unfortunately, the study of the nuclear matter within 
the soliton model often fails to fit the experimental val- 
ues, even in the saturation energy. This may be caused 
by the fact that the topological soliton picture is based on 
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the approximation of large A^c-limit of QCD and therefore 
works weh only in the very low energy scale. Thus our 
model improves slightly the situation in the sense that we 
take into account the quantum correction of 0{1/Nc) to 
the classical soliton mass of 0{Nc). However, it should 
be noted that as our model contains the valence quark 
explicitly, the physical meaning of such Nc counting is 
obscure. Of course, the prescription is still insufficient, 
and the obtained results will still room for improvement. 



II. THE CHIRAL QUARK SOLITON MODEL 

The CQSM was originally derived from the instanton 
liquid model of the QCD vacuum and incorporates the 
non-perturbative feature of the low-energy QCD, sponta- 
neous chiral s ym metry breaking. The vacuum functional 



is defined by 



i ym i 

M 



Z = / VnV^pV'ip^ exp 



where the SU(2) matrix 



(1) 



(2) 



with 



U = e-Kp{iT ■ = —{a + iT ■ tt) (3) 

describes chiral fields, ip is quark fields and M is the dy- 
namical quark mass. We choose the constituent quark 
mass M = 420 MeV which reproduces the experimental 
observables of the free nucleon correctly • /tt is the 
pion decay constant and experimentally /,r ~ 93MeV. 
Since our concern is the tree-level pions and one-loop 
quarks according to the Hartree mean field approach, 
the kinetic term of the pion fields which gives a con- 
tribution to higher loops can be neglected. Due to the 
interaction between the valence quarks and the Dirac sea, 
soliton solutions appear as bound states of quarks in the 
background of self-consistent mean chiral field. Nc va- 
lence quarks fill the each bound state to form a baryon. 
The baryon number is thus identified with the number of 
bound states filled by the valence quarks 2A\ . The B = 1 
soliton solution has been studied in detail at classical and 
quantum level in [ll[ll[ll[lil25|. 

The vacuum functional in Eq. JQ) can be integrated over 
the quark fields to obtain the effective action 

Ses[U] = -iNc\ndet{i^- MW'-) (4) 
= -'-NcSplnD^D (5) 

where iD = — MU""^ {iD is called the Dirac op- 
erator). This determinant is ultraviolet divergent and 



must be regularized. Using the proper-time regulariza- 
tion scheme, we can write 



iNrT 



r ^Spfe-^'^^-e-^iS^"- 

Jl/A2 T ^ 



duj 



dr. 



-Sp 



1/A2 



(6) 



where T is the Euclidean time separation, A is a cut-off 
parameter evaluated by the condition that the derivative 
expansion of Eq. reproduces the pion kinetic term with 
the correct coefficient i.e. 



f 



Nc-m^ 
47r2 



dr 



and H is the Dirac one-quark Hamiltonian defined by 



(7) 



H = 



(8) 



Dq = D(U = 1) and Hq = H{U = 1) correspond to the 
vacuum sectors. At T ^ oo, we have e''^"''' ^ q-^e^^^t ^ 
Integrating over lo in Eq. ^ and constructing a complete 
set of eigenstates of H with 



one obtains the sea quark energy 
Nc f°° dr 



(0) 



(9) 



(10) 



In the Hartree picture, the baryon states are the quarks 
occupying all negative Dirac sea and valence levels. 
Hence, if we define the total soliton energy ^static the 
valence quark energy E^^\^\ should be added; 



(11) 



To obtain the 5 = 1 soliton solution, we impose the 
hedgehog ansatz on the chiral field 

C/(r) — exp{iF{r)f ■ r) = cos F{r) + if ■ t sin F(r) (12) 

with the boundary conditions 

F(0) = -7r, F(cx)) = 0. (13) 

The one-quark hamiltonian (jSJ reads 

H{U'''') = -ia ■ V + /3M(cosF(r) + i-i^f • TsinF(r)) . 

(14) 

This hamiltonian does not commute with the total angu- 
lar momentum j nor the isospin t/2 but commute with 
the grand spin operator K — with K] — 0. H 
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also commutes with the V = "fa which turns to be a par- 
ity operator. As a result, the one-quark eigenstates are 
labeled by the X = 0, 1, 2, • • • and the parity V ^ ±. The 
three valence quarks occupy the lowest states = 0+ 
and are responsible for the baryon number (= 1) (non 
topological charge). 

Field equations for the chiral fields can be obtained by 
demanding that the total energy in Eq. be stationary 
with respect to variation of the profile function F{r), 



5F{r 



■Es- 



which produces 



S{r) sinF(r) = P(r) cos F(r), 



(15) 




free 

R = 2 f m 
R = I fm 
R = 0.5 fm 



2.0 



2.5 



3.0 



where 



S{r) = 7V,^(n,0(e,)-f sign(e,)AA(e,)) 

X (z.|7"<5(|x|-r)|z.), (16) 
P{r) = 7V,^(n,0(e,)+sign(e,)AA(e,)) 



X {i^\ijj''f-TS{\x\—r)\iy) 



(17) 



with 



AA(e,) 



1 ^ (\ /e,A2 



(x) 



and rin is the valence quark occupation number. 

The procedure to obtain the self-consistent solution of 
Ea. (|15|l is that 1) solve the eigenequation in © under an 
assumed initial profile function i^o(^)i 2) use the resul- 
tant eigenfunctions and eigenvalues to calculate S{r) and 
P(r), 3) solve Ea. l|15|l to obtain a new profile function, 
4) repeat 1) — 3) until the self-consistency is attained. 



III. THE NUMERICAL BASIS 

In this section we present the numerical method of 
eigen problem of the hamiltonian l|14|l . The hamiltonian 
with hedgehog ansatz commutes with the parity and the 
grandspin operator given by 

K = j + t/2 = i + (t/2 + r/2, 

where j,Z are respectively total angular momentum and 
orbital angular momentum. Accordingly, the angular ba- 
sis can be written as 



MI^'^3))|iT3) 



FIG. 1: Profile functions for R = 0.5, 1, 2 fm and the free 
soliton. 



For B = 1 solution, following states are possible 



\0) = \iK K + -)KM), 
\1) = \{K K-^)KM), 
\2) = \{K + IK+^)KM), 
\3) = \{K-IK-^)KM). 

With this angular basis, the normalized eigenstates of the 
free hamiltonian in a spherical box with radius R can be 
constructed as follows: 



u 



V 



''KM 

{b) 
KM 

(a) 
KM 

ib) _ 
KM ' 



Nk 
Nk 



iu;XjK{kr)\0) 
^,jK+i{kr)\2) 

ii^tjK{kr)\l) 
-^7jK-i{kr)\3) 

i^tjK+iikr)\2) 
-u;-jK{kr)\0) 

i^tjK-i{kr)\3) 
^7JK{kr)\l) 



(19) 



with 



Nk 



^R'^jK+iikR) 



-1/2 



(20) 



(18) 



and cj+>o,tJ,^<o = sgn(efe), w+ = k/{ek + M). 

The u and w correspond to the "natural" and "unnatural" 
components of the basis which stand for parity (— 1)^ 
and (—1)^"*"^ respectively. The momenta are discretized 
by the boundary conditions jK{kiR) = 0. The orthogo- 
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TABLE L The classical mass for the original Kahana-Ripka 
basis and modified version (in MeV), with M — 400 MeV, 
R — 6 fm. The error becomes of order ~ 10"''. 





Valence 


Vacuum 


Total 


free 


191 


637 


1209 


modified 


192 


633 


1210 



nality of the basis is then satisfied by 



drr'^ jK±i{kir)jK±i{kjr) 



(21) 



Let us examine the boundary conditions for the chiral 
and Dirac fields to construct the nuclear matter solu- 
tion in the Wigner-Seitz approximation. When the back- 
ground chiral fields are periodic with lattice vector a, the 
quark fields would be replaced by Bloch wave functions 
as il){r + a) = e^^'°"il){r). In the Wigner-Seitz approxi- 
mation, however, the soliton is put on the center of the 
spherical unit cell with the radius R [a ~ 2R) and the 
dispersion k is assumed to be zero. For the profile func- 
tion F{r), the periodicity and the unit topological charge 
inside the cell require the boundary conditions 



a'(0) 
7r(0) = 



-a'{R) 



= 




F(0) = -7r,F(i?) =0. (22) 



For the Dirac eigenstates, modification in the basis is 
needed. For odd number of K, the boundary condition 
is same as the free case with 

jKikiR)^0. (23) 

For even K, the following conditions must be satisfied 



jK+i(fci"^i?) 
JK-iikl'^R) 



0, 
= 0, 



for 
for 



(a) (a) 

(b) (b) 
"^KM^ ""km ■ 



(24) 



Obviously the conditions (|24l) partially break the orthog- 
onality of the basis (|21|l for the finite value of R. However 
we can solve the eigenvalue problem properly (see Ta- 
ble QJ. Although the motivation is different, the similar 
treatment has been already introduced in Ref . 26] . 

Fig. 1 shows the self-consistent profile functions for 
free (i? oo) and various values of the cell radius R. 
In Fig. 2, we present the results of the classical energy 
of the soliton and its valence and vacuum contributions 
as functions of R. We find the shallow minimum of the 
classical energy at R ^ 1.2 fm. 



IV. SPURIOUS CENTER OF MASS 
CORRECTION 
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FIG. 2: Classical soliton energy and its valence and vacuum 
contributions lllll . 



The minimum found in Fig. [21 is not regarded as a true 
saturation point because it contains the zero-point en- 
ergy contributions. The quark contribution to the mean- 
field expectation value of the square of the total momen- 
tum appears at the classical level although it should 
be zero because the soliton is rest at the cell center in 
the present approximation. Therefore, the correspond- 
ing kinetic energy should be subtracted from the total 
energy. The effects of the spurious center of mass mo- 
tion is roughly estimated by the method of Pobylitsa et 
al. [231 . Considering the translational degrees of free- 
dom and performing their quantization, one obtains the 
correction at a rest frame as 



static 



'static 



= E. 



'static 



(25) 



The correction is easily evaluated by using the numerical 



basis given in Eq. H19(l as P 



.2 (o) 

"■KM 



{hr) 



(a) 
''KM 



ihr). 



As is shown in Fig. |3 the minimum disappears after re- 
moving the zero-point energy contributions (|25|l . This 
is explained by the observation that the contribution of 
the center of mass motion becomes small with increasing 
density (see Fig. ^ and the caption) . 



V. COLLECTIVE QUANTIZATION 

The solitons that we have obtained in the previous 
section are classical objects and therefore must be quan- 
tized to assign definite spin and isospin to them. For the 
solitons in the free space, quantization can be performed 
semiclassically for their rotational zero modes. For the 
hedgehog sohton, because of its topological structure, a 
rotation in isospin space is followed by a simultaneous 
spatial rotation. Let us introduce the dynamically ro- 
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R[im] r[fm] 



FIG. 3: Classical soliton energy after removing the spurious 
center of mass motion 1251 . 



tated chiral fields 



U{x, t) = A{t)U{x)A{t)\ A{t) e SU(2) 



(26) 



In a crystal configuration, the solitons are fixed on the 
spatial lattice point and their isospin orientation is chosen 
so as to minimize the eirergy of the system. If one rotates 
each soliton in the crystal, it changes the isospin orien- 
tation and increases the energy. Thus there is only one 
isospin collective coordinate corresponding to the over- 
all orientation of the crystal in isospace, called global 
isospin, in the soliton crystal [lll2^. 

The Wigner-Seitz treatment with spherical cell approx- 
imation may cure the situation. Because in this ap- 
proximation the information of the crystalline structure, 
hence, the isospin structure is completely lost at least 
in the low-density, the rotational zero-mode would be 
recovered. Thus, we apply the zero- mode quantization 
method to the WS-cell to estimate the nucleon and the 
delta mass splitting in the matter. 

By transforming the rotating frame of reference, the 
Dirac operator with Eq. (|26|l can be written as 



where 



ib = A{t)-f°[idt - H{U'"') + n]A{t)^ 



2 



(27) 



(28) 



n is the angular velocity operators for an isorotation. 
Assuming that the rotation of the soliton is adiabatic, we 
shall expand the effective action SeS around the classical 
solution U{x) with respect to the angular momentum 
velocity SI up to second order 



-iN.Sp [\og{idt' 



H + n)] -Sp[logiidt~H)] 



FIG. 4: The "upper" u(r) and the "lower" w{r) component 
of valence quark wave functions for various cell radius R with 
the boundary condition 'w{R) = 0. Non vanishing values of 
upper component at the cell boundary u{R) come from the 
zero-mode elements in the basis. 



With the proper-time regularization, we have 

ah 

where /sea. ah is the vacuum sea contribution to the mo- 
ments of inertia defined by 

/soa,ah = ^N^^f{€^,e^,^.){v\Ta\^i){^i\n\v) , (29) 

with the cutoff function /(e^, e^^ A) 

2Ae2^^^^^e2^ 

^ sgn(e^)erfc(|e^|/A) - sgn(et,)erfc(|et.|/A) 

Similarly, for the valence quark contribution we have 

(val|Ta|/i) (/i|rh| val) 



(30) 



/^^val 



The total moments of inertia are then given by the sum 
of the vacuum and valence, lab = Iva.i,ab + /sea,ab- The 
hedgehog ansatz of the chiral fields ensure the relation 
for the moment of inertia 



hi — I22 — I: 



33- 



(31) 



The quantization condition for the collective coordinate, 
A{t), define a body-fixed isospin operator K as 

I.^n'^-trfA^-^^^-ka, (32) 
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FIG. 5: Moment of inertia : the vacuum H29|l and the valence 
1301 1 contribution and their sum. 



FIG. 6: Quantized soliton energies of nucleon A'^ and delta 
resonance A(1232) lIMIl . 



These are related to the usual coordinate-fixed isospin 
operator ia by transformation, 



1 



TT[TaA{t)T''A{t)^]h 



(33) 



To estimate the quantum energy corrections, let us in- 
troduce the basis functions of the spin and isospin oper- 
ators which were inspired from the cranking method for 
nuclei 



{Aliisks) 



2i + l 
87r2 



i-iy+^^D'(A) 



where D is the Wigner rotation matrix. Finally, we find 
the quantized energies of the soliton as 



E — i?si 



21: 



33 



(34) 



where i{i + 1) is eigenvalues of the Casimir operator i^. 
The moment of inertia for the vacuum H29I) and valence 
(|30|l and their sum are given in Fig. O In Fig. we 
present the energy of nucleon (i = i) and A (i = |). 

In this cranking procedure, the zero-point energy of 
the rotational motion (i^)/2/33 must be removed from 
Eq. 131 53,E3- Finally, we obtain the mass of nucleon 
and delta 



En — -E'static " 
Ea = -2'static 



4/33 

3 



4/: 



33 



(35) 
(36) 



Fig. 13 shows the energy of nucleon and delta after sub- 
tracting the spurious zero-point energy. The minimum 
for nucleon is observed at i? ^ 1.8 fm which corresponds 
to the density ~ 0.04 fm^'^. This value is much 



lower than the experimental value. The binding energy 
is Eb ~ 18 MeV which is not far from the experimen- 
tal observation. For A, we also find the shallow minima 
at i? ^ 1.22 fm which corresponds to pA ~ 0.13 fm^'^. 
The A saturation is attained at the density pa/pn ~ 3.2 
which is close to the prediction of density pa/pn ~ 2 — 3 
in the framework of the quantum hadrodynamics 32, 331. 
The advantage of our approach is that the model does 
not require any tuning parameter for the A spectra in 
the hadrodynamics calculations. 



VI. SUMMARY 

We have studied soliton solutions in nuclear medium by 
using the Wigner-Seitz approximation. The chiral quark 
soliton model was used to obtain the classical soliton so- 
lution. In this letter we especially focused on the prop- 
erties of nucleon and A in matter. We quantized the 
soliton semiclassically. The adiabatic rotation for the 
(iso-)rotational zero mode was performed and the nu- 
clear saturation points were obtained for nucleon and A 
matter. 

Here we did not consider the following effects which 
should be investigated to develop our understanding of 
the dense nuclear matter: 

• band structure of the quarks 

• R dependence of the constituent quark mass M and 
the cutoff parameter for the vacuum A 



• inclusion of the heavier mesons {p,uj, 
soliton solutions 



to the 



improvement of the correction by the zero point 
energy and Casimir effects 
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FIG. 7: Masses of A'^, A, after spurious energy subtrac- 
tions 135l l- ll36l l. 
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FIG. 8: Binding energy of nucleon for the various constituent 
quark mass M (in MeV). 

• quark-meson couplings and the Fermi motion of the 
baryons 

• crystaUine order in high density phase 

• SU{3) extension. 

As is expected, our model provides much lower value of 
the saturation density than the experiment. In this anal- 
ysis, the Wigner-Seitz cell is approximated by a sphere 
and thus high density matter is attained by shrinking the 
cell volume with the spherical shape of each soliton un- 
changed. However, in reality, the neighborhood solitons 
start to overlap and the structure will be deformed from 
uniform nuclear matter at high density. In this phase, the 
hedgehog ansatz should not be appropriate any more. 



FIG. 9: Kinetic energy of the pion 



We observed the increase in the zero mode of the cen- 
ter of mass motion of the soliton for higher density, which 
means that the soliton tends to rest in the WS approxi- 
mation. In this case, we should employ the exact WS cell 
which reflects the background crystal symmetry instead 
of sphere, to get higher saturation density. The inclusion 
of band effects may also improve our results. In Ref. W^, 
the authors imposed the Bloch-like boundary conditions 
on the s-wave valence quark wave function and estimated 
the soliton energy self-consistently. They found that the 
effects of the admixtures of higher states are small except 
for the scalar quark density. In fact, the band structure 
will appear at some critical density and the correction for 
the quantum energy may become more important at the 
dense medium because the radius of the soliton, that is, 
the moment of inertia, strongly depend on the position 
of the band l^. 

Generally speaking, the constitue nt q uark mass M is 
momentum- and density-dependent [l6j- We chose the 
value M = 420 MeV as it reproduces the free nucleon 
observable. In Fig.|Sl one can see that for larger value of 
M, the saturation point goes to inward and the binding 
becomes deeper. Varying the value of M for each density 
may give a better result for the saturation point. 

An important feature of the nucleon in a matter is 
about its size. It is believed that the nucleon will swell in 
the medium. The authors of Ref. [s^ observe such effect 
with reducing effective quark mass M* in the Nambu- 
Jona-Lasinio type quark-soliton model. We confirmed 
within our model that as smaller the M, the size of the 
soliton increases. But in that case, the saturation be- 
comes shallow (Fig.|Sll. Recently, we investigated soliton 
solutions in the CQSM taking into account p, w mesons 
which will improve the short distance physics. We are 
able to obtain deeper binding energy as decreasing the 
value of M. We will report it on forthcoming article. 

In Fig. El one finds that the spectra of nucleon and 
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A are too small compared to the experimental values. 
Obviously it is due to the subtraction of the zero-point 
corrections. A little more sophisticated approach of the 
spurious motion is performed in Ref. and by applying 
this approach to our analysis, the results will be improved 
to a certain extent. Also, the meson coupling to the 
quark inside nucleon and A should be important to shift 
the minima at higher density. 

In Fig. Eland Fig.[7| one finds the nucleon- A mass dif- 
ference gradually decreases as matter density increases 
and eventually it vanishes. The reduction in the mass 
difference has been observed previously in a similar chi- 
ral soliton model but employing somewhat different pro- 
jection technique for quantum number [^MUBl- In the 
present formulation, the behavior is not fully under- 
stood because it should be explained by the dynamics 
of hadrons, that is, QCD. In the naive SU{6) quark 
model, the mass difference is ascribed to the hyperfine 
splitting 37]. The reduction may imply the increase of 
the distance between quarks. In fact, in Fig. ^ one can 
see the concentration of the quarks at the cell boundary 
as the density increases. 

Alternatively, if we understand the A as a compos- 
ite object (resonance state) of the nucleon and pion, the 
mass difference can be interpreted as the energy of pi- 
ous bound to the nucleon. Although it is absent in the 



present formulation, the pion kinetic energy inside the 
soliton can be estimated as 

In Ref. |22| , the authors introduced the r- and the cutoff 
parameter of the vacuum A-dependent form of the pion 
decay constant /^r (r. A) and estimated its density depen- 
dence with the A whose value is set for the free space 
value of /tt . The (r, A) determined in such a way is 
essentially valid only for the free space limit R oo. 
Therefore we shall simply take the value in free space 
/tt — 93 MeV. Fig. shows the result of the kinetic en- 
ergy of pions as a function of R and one can observe that 
the energy is reduced as the density increases. This re- 
duction of the pion kinetic energy may contribute to the 
reduction of the mass difference. 

Our formulation is directly applicable to the SU (3) 
octet-decuplet baryon spectra in nuclear matter l38j | . 
After the above effects are properly incorporated and 
more realistic estimation of the saturation points is 
achieved, it will be interesting to study the SU{3). 



[1] Igor Klebanov, Nucl. Phys. B 262, 133 (1985). 

[2] E. Wiist, B. E. Brown and A. D. Jackson, Nucl. Phys. A 

468, 450 (1987). 
[3] Alfred S. Goldhaber and N. S. Manton, Phys. Lett. B 19, 

231 (1987). 

[4] L. Castillejo, P. S. Jones, A. D. Jackson, J. J. M. Ver- 
baarschot and A. Jackson, Nucl. Phys. A 501, 450 (1987). 

[5] M. Kugler and S. Shtrikman, Phys. Rev. D 40, 3421 
(1989). 

[6] Byung-Yoon Park, Dong-Pil Min, Mannque Rho and 

Vincente Vento, Nucl. Phys. A 707, 381 (2002). 
[7] Joachim Achtzehnter, Werner Scheid and Lawrence 

Wilets, Phys. Rev. D 32, 2414 (1985). 
[8] B. Banerjee, N. K. Glendenning and V. Soni, Phys. Lett. 

B 155, 213 (1985). 
[9] N. K. Glendenning and B. Banerjee, Phys. Rev. C 34, 

1072 (1986). 

[10] Detlev Hahn and Norman K. Glendenning, Phys. Rev. C 

36, 1181 (1987). 
[11] Urban Weber and Judith A. McGovern, Phys. Rev. C 

57, 3376 (1998). 
[12] H. Reinhardt, B. V. Dang, and H. Schuiz, Phys. Lett. B 

159, 161 (1985). 
[13] M. C. Birse, J. J. Rehr and L. Wilets, Phys. Rev. C 38, 

359 (1988). 

[14] Nir Barnea, Timothy S. Walhout, Nucl. Phys. A 677, 
367 (2000). 

[15] M. Kutschera, C. J. Pethick and D. G. RavenhaU, Phys. 

Rev. Lett. 53, 1041 (1984). 
[16] D. I. Diakonov, V. Yu. Petrov, and P. V. Pobylitsa, Nucl. 

Phys. B 306, 809 (1988). 



[17] H. Reinhardt and R. Wiinsch , Phys. Lett. B 215, 577 
(1988). 

[18] Th. Meissner, F. Griimmer, and K. Goeke, Phys. Lett. B 

227, 296 (1989). 
[19] For detailed reviews of the model see: 

R. Alkofer, H. Reinhardt and H. Weigel, Phys. Rept. 265, 

139 (1996); 

Chr. V. Christov, A. Blotz, H.-C.Kim, P. Pobylitsa, T. 

Watabe, Th. Meissner, E. Ruiz Arriola, K. Goeke, Prog. 

Part. Nucl. Phys. 37, 91 (1996). 
[20] D. Diakonov, V. Petrov and M. Polyakov, Z. Phys. A 

359, 305 (1997). 
[21] T. Nakano et al, Phys. Rev. Lett. 91, 012002 (2003). 
[22] P. Amore and A. De Pace, Phys. Rev. C 61, 055201 

(2000). 

[23] I. Adjah, I. J. Aitchison, and J. A. Zuk, Nucl. Phys. A 

537, 457 (1992). 
[24] S. Kahana and G. Ripka, Nucl. Phys. A 429, 462 (1984). 
[25] M. Wakamatsu and H. Yoshiki, Nucl. Phys. A 524, 561 

(1991). 

[26] H. Weigel, R. Alkofer and H. Reinhardt, Nucl. Phys. B 

387, 638 (1992). 
[27] P. V. Pobylitsa, E. Ruiz Arriola, Th. Meissner, F. 

Griimmer, K. Goeke and W. Broniowski, J. Phys. G 18, 

1455 (1992). 

[28] W.K.Baskerville, Phys. Lett. B 380, 106 (1996). 

[29] L. C. Biedenharn, Y. Dothan and M. Tarlini, Phys. Rev. 

D 31, 649 (1985). 
[30] Thomas D. Cohen and Wojciech Broniowski, Phys. Rev. 

D 34, 3472 (1986). 
[31] A. Bohr and B. Mottelson, Nuclear structure, Vol.11 



9 



(World Scientific Publishing Co. Pte. Ltd, Singapore, 

1998). 

[32] B. M. Waldhauscr, J. Tliois, J. A. Maruhn, H. Stocker, 
and W. Greiner, Phys. Rev. C 36, 1019 (1987). 

[33] Zhuxia Li, Guamgjun Mao, Yizhong Zhuo, and Walter 
Greiner, Phys. Rev. C 56, 1570 (1997). 

[34] Chr.V.Christov, K.Goeke, Nucl. Phys. A 564, 551 
(1993). 

[35] E.Ruiz Arriola, Chr.V.Christov and K. Goeke, Phys. 



Lett. B 225, 22 (1989). 
[36] Chr.V.Christov, M. Fiolhais, E.Ruiz Arriola and K. 

Cocke, Phys. Lett. B 243, 333 (1990). 
[37] A. De Rujula, H. Georgi and S. L. Glashow, Phys. Rev. 

D 75,147 (1975). 
[38] A. Blotz, D. Diakonov, K. Goeke, N. W. Park, V. Petrov 

and P. V. Pobylitsa, Nucl. Phys. A 555, 765 (1993). 



